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ABSTRACT 

In this Letter, we give new constraints on planet migration. They were obtained under the assumption 
that Saturn's current obliquity is due to a capture in resonance with Neptune's ascending node. If planet 
migration is too fast, then Saturn crosses the resonance without being captured and it keeps a small 
obliquity. This scenario thus gives a lower limit on the migration time scale r. We found that this 
boundary depends strongly on Neptune's initial inclination. For two different migration types, wc found 
that r should be at least greater than 7 Myr. This limit increases rapidly as Neptune's initial inclination 
decreases from 10 to 1 degree. We also give an algorithm to know if Saturn can be tilted for any migration 
law. 
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Hl^ Introduction 

^ ■ It is now well accepted that the Solar System was 
I Tin 'orn compact after the protoplanetary gas disk dis- 
sipated. Then planets migrated due to interactions 
^~ ^ith the primordial planetesimal belt. The Nice model 
>TlGomes et al.ll2005l: iTsiganis et al]|2005l: iMorbideUi et all 
gives a unified scenario of this planetary migra- 
^^ibn, but it is still not fully constrained. For example, 
^— ihe Nice model allowed two possible classes of late evolu- 
^T;ion (fNcsvorny et al. 2 00^ . In the first one, called "class 
0\lA" , Neptune is scattered to 22-25 AU and reaches its 
orbit by slowly migrating over more than 5 AU. In 
^^he second class, labeled "DE", Neptune is placed to its 
TTiirrent orbital distance with large eccentricity w 0.3 and 
J>then slowly circularizes. Besides, the Nice model does not 
Il^^cJnstrain inclinations, and the time scale of this late evo- 
>^Ution is uncertain. Nevertheless, constraints on the mi- 
j^ration time scale were obtained from the distribution of 
■ " the K uiper belt on the one hand ()Murrav-Clav fc Chianal 
[20051) . and fr om the distribution of the main asteroid belt 
on the other (jMinton fc Malhotral [20091 ) . They both as- 
sumed an MA migration type without long-term evolution 
of eccentricities and inclinations. The former obtained a 
migration time scale r between 1 and 10 Myr, and the 
latter found r < 0.5 Myr. 

In this Letter, we aim at giving new constraints based 
on Sa turn tilting (jWard fc IIamiltonll2004HHamilton fc Ward! 
|2004[ ). According to Ward and Hamilton, Sat urn's large 
obliquity, e = 26.73919 deg ([Helled et al.ll2009t) . is due to 
a resonance capture between its spin axis and Neptune's 
orbit. Given the large uncerta inties on Saturn's preces - 
sion rate -0.75 ± 0.21"- yr"! ([Ward fc Hamiltonll200l . 
the more accurate regression of Nep tune's orbit plane 
sg = -0.692"- vr" ^ (iLaskar et al.ll2004f) is in deed included 
in the errorbars. [Ward fc HamiltonI (|2004l ) assume that 



today the two frequencies are equal. In their scenario, 
the norm of the frequency of Neptune's ascending node 
was initially larger, and then it captured Saturn's spin 
axis as it decreased due to Neptune's migration and/or 
the d issipation of the planetesimal disk ([Ward fc Canud 
l2006t ). In their numerical model, they took a quasiperi- 
odic model of the Solar System and forced an exponential 
evolution of the frequency sg- Here, we show that Saturn 
can tilt in both migration classes, and that it gives a lower 
limit on the migration time scale. This limit depends on 
Neptune's initial inclination. 

A recent paper bv lHelled et"an(|2009D seems to contra- 
dict Ward and Hamilton's scenario. It gives a new esti- 
mate of Saturn's precession rate — 0.7542 ± 0.0002"- yr~^ 
that is incompatible with a resonance with sg. We show 
that with this value, Saturn can still evolve to its current 
state but that it is very unlikely. We discuss this result 
in our conclusion. 

2. Spin axis evolution 

Here we recall the equations of motion of a planet axis 
and give the current dynamical state of Saturn's spin axis. 
The evolution of the spin axis w of a planet in a fixed 
reference frame (i,j,k) where k is the direction of the 
total orbital angular momentum, is given by 



dw 
~dt 



= — Q'(n ' 



w n X w 



(1) 



where n = *{nx, ny,nz) is the normal to the orbit, and a 
is the precession constant. Without planetary perturba- 
tions, n is fixed and the spin axis w precesses uniformly 
around n with constant obliquity cos e = n - w. How- 
ever, in a multi-planetary system, n evolves due to sec- 
ular interactions. The long-term evolution of n can be 
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approximated by a quasiperiodic expression 



(2) 



where the Vk (sorted with increasing amphtudes Ik) 
are corabinat ions of the fundamental frequencies gj, s,- 
(|Laskaij[T99l . For Saturn, 1^2 = ss = -0.692" • yr""^ 
and I2 = 0.064 deg. As the other terms have only 
very weak effects on the behavior of Saturn's spin axis 
([Hamilton fc Wardll2004D . one can retain this single term 
in the orbital precession, which makes the problem inte- 
grable (|Colombdll966H H enrard fc Murigand3ll987t ). The 
associated autonomous Hamiltonian, written in a moving 
reference frame related to the orbital plane, reads 



(3) 



with now k ^ *(/,0, Vl - P) and n = *(0,0,1). The 
equation of motion, obtained from dm/dt = x w, 

is 

r/w 

(4) 



= —a(n ■ w)n x w — z/k x w . 

dt 

This system possesses 4 relative equilibriums named 
Cassini states for which the three axes w, n and k are 
collinear, and a separatrix delineating 3 zones in the phase 
space (see Fig.[T]). Hereafter, we label the 3 zones after the 
Cassini state they contain. Saturn's spin axis coordinates 
in the orbital frame are 



^sin e cos V'^ 
sin e sin ip 
< cos e J 



(5) 



with e = 26.73919 deg (iHelled et all I2OO90 and ip = 
-31 deg ([Hamilton fc WardI I2004D . As ip ^ 0, the sys- 
tem is not in a Cassini state. Given th e coordinates of w, 
if -a cos e G [-0.730, -0.666] " • yr'^ (|Ward fc HamiltonI 
I2004D then Saturn's spin axis is in resonance around 
Cassini state 2 with a libration amplitude larger than 
31 deg (see FigHb), else it is in circulation around either 
Cassini state 1 (see FiglT^) or Cassini state 3. Literature 
gives three different values of S aturn's p recession rate. 
Two of them, -0.74± 0.7"- vr"^ (iFrench et al 19931) and 



-0.75 ± 0.21"- yr"^ (|Ward fc HamiltonI |2003) are com- 
patible with a libration in zone 2, whereas th e third one, 
-0.7542 ± 0.0002"- yr-i (IHelled et all [20091) . constrains 
Saturn's axis to circulate in zone 1. In the following, we 
study these two c ases. In Cas e I, we use the precession 
constant given bv IHelled et al.l (j2009f ). and in Case II we 
set a such that —a cose = = -0.692"- yr~^. In our 
numerical integrations detailed below, we take into ac- 
count the dependences of a in Saturn's semi-major axis 
and eccentricity. 

3. Orbital evolution 

We integrate the sec ular equations of motion derived 
from the Hamiltonian of lLaskar fc Robutell (|1995l ) written 
up to degree 4 in inclinations and eccentricities. In or der 
to fit to the present value of ss (iLaskar et al.ll200l ). a 
small constant offset 6ss = —0.00342"- yr~^ is added in 
the model. This offset was obtained by frequency analysis 




Fig. 1. — Projection of the spin axis v^f ([5|) in the orbital 
frame (in abscissae sinecos'i/', in ordinate sin e sin -0 ) . (a) 
Case I with a = 0.845"- yr"! (IHelled et all 12009 ). (b) 
Case II with a = 0.775"- yr'^ (|Ward fc HamiltonI I2004D . 
Cassini state 3 corresponds to a retrograde rotation of 
Saturn and is not represented in these figures. The current 
position of Saturn's spin axis is represented by a large 
filled circle. The small filled circles are Cassini states and 
the curves are energy contours. The bold curve is the 
separatrix that delineates the libration area in grey. 



(|Laskarlll990f ) of our analytical model (Table [T|). For the 
class "MA", we consider only the last 3 AU migration 
of Neptune. When Neptune was closer to the Sun, the 
frequency sg was too large to have any effect on Saturn's 
axis. Migration is simulated by an additional force leading 
to the following exponential law, 



ait) = ao + Aa(l - e"*/^), 



(6) 



with Aa = +0.1,-0.3,-1.3,-3 AU respectively for 
Jupiter, Saturn, Uranus, a nd Neptune. It is scaled from 
iMinton fc Malhotra (l2009l) and i t is in a greement with the 
full integration of lTsiganis et al.l (|2005f ) . In the same way, 
for the class "DE" we apply an external force that gives a 
long-term exponential evolution of Neptune eccentricity 
starting at cq = 0.3 and finishing at its current value. 
For both classes, we did integrations with constant Nep- 
tune inclination, and others with an exponential damping 
with the same r. For each value of r, an integration in 
the past is done to obtain initial conditions for the or- 
bital coordinates. Saturn's initial obliquity is then set to 
Co = 1.5 deg. 



Table 1: Secular Frequencies Associated with the Preces- 
sion of the Ascending Nodes. 



Frequency 


Laskar et al. (2004) 


Secular integration 


S5 


-0.000 


-0.000 


S6 


-26.348 


-26.569 


S7 


-2.993 


-2.996 


S8 


-0.692 


-0.689 



We now look at the effect of the dissipation of the re- 
maini ng primordial plan etesimal belt. Following a sugges- 
tion of lMorbidellil (|2009l ). the mass rriK of the planetesimal 
belt in the class "MA" is estimated by energy conservation 
as f ollows. Initially, planet esimals are distributed follow- 
ing iMorbidelli et al.l (|2004l . fig. 1) and during Neptune 
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Fig. 2. — Results for the MA type migration. Projec- 
tion of Saturn's spin axis on the invariant plane in a 
frame rotating at Neptune regression frequency (in ab- 
scisse sin6'cos(</) — f2), in ordinate sm9sm{(j) — fl)) (see 
Eq. [9]). The filled circle represents its current position. 
Subfig. la, lb and Ic. Case I with r = 180 Myr and 
1/; = 108.674, 108.675 469, 108.675 475 deg. Subfig. Ila, 
lib and lie. Case II with r = 20,200,300 Myr and 
ip = deg. 



migration, planetesimals move from their initial position 
to Uranus' orbit. This leads to niK = 1.7 ± O.IM®. Dur- 
ing planet migration with a planetesimal disk, we force an 
exponential decrease of the planetesimal belt mass with 
the same time scale as the semi-major axis one. To first 
order, the averaged effect of a planetesimal of mass rrii and 
semi-major axis on Neptune's nodal precession rate is 

(5s8(mj,ai) = b''y\{aN/ai) (7) 

with njv and being respectively Neptune mean motion 
and semi-major axis, and toq the mass of the Sun. We 
model the planetesimal belt by a single annulus with semi- 
major axis aK such that 

Sssinik,ak) = ^ Sssinii, a^) . (8) 
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Fig. 3. — Probability that Saturn librates in zone 2 with 
an amplitude larger than 31 deg as a function of the mi- 
gration time scale r in Case II. MA migration type in grey 
and DE migration type in black. 



Using the mass distribution of iMorbidelli et al.l ()2004[ ) , 
we found a/; = 60 ± 5AU. We run numerical integrations 
with and without a 2 Earth mass primordial planetesimal 
belt and found that the constraints on the migration time 
scales were unchanged. In the following, we give only the 
results of our integrations without a planetesimal disk. 

4. Results 

In Case I, Saturn's spin axis circulates around Cassini 
state 1 with a large obliqu ity (Fig. [1^). To show whether 
it is compatible with the IWard fc Hamilto n (2004) sce- 
nario or not, we did several numerical integrations with- 
out long-term evolution of Neptune inclination and enu- 
merated those ending in zone 1 with an obliquity larger 
than or equal to 26.73919 deg. We first considered the 
class MA of the Nice model and varied the migration 
time scale r from 100 Myr to 600 Myr every 10 Myr. 
Then, for each value of r, we searched the range of the 
initial precession angle t/j for which the final state cor- 
responds to our criterion. We found only 3 values for 
T satisfying the criterion: r e {170,180,190} Myr. In 
each case, the range of possible values for ijj is extremely 
small < 10~^ deg (Fig. [21 labc). Thus assuming an 
equiprobable initial phase, the probability to find Saturn 
in its current state through this mechanism is less than 
3 X 10~^ for any of the three selected r. With a DE mi- 
gration type, the widths of the initial longitude intervals 
are identical. The only changes are in the values of 
the migration time scale r leading to the large obliquity 
circulation state: r e {150, 260, 290, 310, 320} Myr. We 
discuss the implications of these results in the conclusion. 

In Case II, Saturn spin axis is presently in resonance 
with Neptune's ascending node. In that case, planet 
migration must be slow enough for the capture to oc- 
cur, but if it is too slow, then the evolution becomes 
adiabatic and the libration amplitude is too small (less 
than 31 deg) (Fig [21 He). This latter constraint disap- 
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Fig. 4. — Minimal migration time scale as a function 
of Neptune's initial inclination for both migration types: 
MA (solid line) and DE (dashed line). Bold curves are 
results of numerical integrations. Thin curves were ob- 
tained by the algorithm described at the end of Section [5l 



S )ears if the pre capture obliquity is larger than 4.5 deg 
Ward fc Hamilton ,20041 ). We performed 2100 integra- 
tions for each of the two migration types MA and DE, r 
going from 10 to 600 Myr every 10 Myr and ip between 
and 350 deg every 10 deg. The results are summarized 
in Fig. [3l MA type in grey and DE type in black. Prob- 
abilities are now significant and reach 1 for a few time 
scales. We see clear lower limits, r > 90 Myr (resp. 
T > 170 Myr) for the MA (resp. DE) migration type. 
The difference in the results between the two Nice model 
classes comes mainly from the different dependence of the 
semi-major axis and the eccentricity on Neptune's regres- 
sion frequency. In all these integrations, Neptune's incli- 
nation docs not undergo long-term evolution. However, 
the amplitude I2 of Saturn orbital quasiperiodic motion 
([2) is proportional to Neptune inclination. In Section [51 
we show that the higher the inclination amplitude is, the 
faster a planet can be tilted. We thus studied the min- 
imum time scale, for which Saturn's axis ends in zone 2 
with a libration amplitude larger than 31 deg, as a func- 
tion of Neptune's initial inclination (Fig. SI bold curves). 
In both migration classes, Tmin decreases rapidly to ~ 20 
Myr when Neptune's initial inclination increases to 4 deg 
and then it decreases slowly down to w 7 Myr when Nep- 
tune's inclination goes to 10 deg. 

5. Fastest tilting 

In this section, we compute analytically the minimal 
time required to tilt a planet as a function of its inclination 
I{t). We give also an algorithm to check whether Saturn 
can be tilted or not for a given migration. 

We call 9 the inclination of a planet equator relative 
to the invariant plane. As Saturn current inclination is 
small relative to its obliquity e, the two angles 9 and e are 



similar. Let $ and fl be the longitude of the ascending 
node of the equator and of the orbit in the invariant plane. 
We have 



sin 6* sin $ ^ 
- sin 9 cos $ 
cos 9 / 



(ij,k) 



' sin / sin ^ 
— sin / cos n 
. cos / / 



(ij.k) 



and 



sin 9 sin / cos(<I' — fi) + cos 9 cos /, 



(n X w) • k = — sin6'sin/sin(<l> — Q). 
Thus, from (P), 
dcos ( 



(9) 



(10) 



dt 



a sin 9 cos 9 sin / cos /[I +7 cos(<I' — il)] sin($ — il) 

where 7 = tan 9 tan / can vary from to infinity depend- 
ing on the value of the obliquity. We now choose H. that 
maximizes this time derivative as a function of / and 9. 
Doing so, we ensure that it is not possible to have a faster 
evolution of the equator inclination 9. This leads to 



cos($ - il) 



-1 + v/l + 87^ 

47 



(12) 



Substituting this expression in (jlip gives the maximal 
speed 6max such that d9/dt < 8max 



a 



V2 
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cos sin 2/1 



(v/TT8^ + 3)'^' 

7 , \ 1/2 ■ 

(v/TTV + ij 



(13) 



After some calculus, it can be shown that 0max is an 
increasing function of tan/. Thus, if the only constraint 
on the orbit inclination amplitude is an upper limit /max < 
7r/2, the fastest evolution is obtained for / ~ /max- In two 
asymptotic cases, the expressions of 0max arc simpler. For 
/ < \tt/2 - 6*1 or < \tt/2 - /|, we have I7I < 1 and thus 



I cos 6'sin 2/1 



(14) 



In the other case, if \9 - 7r/2| < / or |/ - tt/2\ < 9, the 
parameter 7 is arbitrarily large and (|13p becomes 



— sin 6* sin /. 
2 



(15) 



Using the approximation for small angles I7I ^ 1 Hi 
the minimum time imin required to bring 9 from to ^ond 
at constant inclination amplitude / is 



a sin 21 



■In 



tan 



^cnd 



(16) 



In Saturn's case, the amplitude of the mode responsible 
for the tilt is /2. Whenever Neptune's inclination is less 
than 10 deg, I2 remains below 0.9 deg, 7 < 8 x 10"'^ and 
Expressions [TB)) are valid. This minimum time tmin 
decreases with the inclination amplitude /. For example, 
with I = I2 ~ 0.064 deg. Equation ^6)) gives tmin = 105 
Myr in Case I and tmin = 115 Myr in Case 2. 
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From this study, it is possible to check whether Sat- 
urn's axis can be tilted or not for a given migration. Let 
6*2 be the value of 6 at the Cassini state 2. O2 and its 
time derivative O2 = 82 are functions of orbital parame- 
ters through /(t), v{t) = ss{t) and a{t). During a tilt, 9 
oscillates around the increasing 62 and = 6 reads 



6 = 62 + Aujiib sin(u)iib< + 



(17) 



where A and ip are respectively the libration amplitude 
and a phase, and a ; iib is the libration amplitude given by 
iHamilton fc WardI (|2004D 



wiib 



\/— ai^sin^sin/. 



(18) 



For a given migration, Saturn's axis can tilt if and only if 
there exist A and <f> such that Q < Omax (|14l and[T7 |l dur- 
ing all the evolution. Replacing 9 by 62 in p^andfTSjl. one 
obtains a criterion that depends only on orbital parame- 
ters. We applied this criterion on the systems studied in 
Section m For each value of Neptune's initial inclination, 
we integrated once the system with a given r. Then, we 
rescaled the derivatives 9 for different value of r until 
the criterion is verified. The resulting values of T,„i,i are 
displayed in Fig. (jH thin curves). 

6. Conclusions 

First of all, we see that the lHelled et"an (1200911 preces- 



sion c onstant is incompatible with the lWard fc Hamilton 
(2004) scenario. This is a robust result. iHelled et al 



20091 ) obtained Saturn's precession constant from an em- 
pirical model of its internal structure. They used Saturn 
mass, radius and gravitational coefficients J2, J4, and Jg 
to fit a density profile represented by a sixth degree poly- 
nomial. From this density profile they derived the normal- 
ized axial moment of inertia 7 directly related to the pre- 
cession constant. Our results suggest, rather, considering 
7 as an additional independent parameter to better con- 
strain Saturn's interior. If Saturn is a ctually in libration 
in zone 2 then 0.2257 < 7 < 0.2438 (|Ward fc HamiltonI 

I2004D . 

Assuming the lHamilton fc WardI (|2004[ ) precession con- 
stant, Saturn's spin axis is likely to evolve toward a libra- 
tion in zone 2 whatever the migration class is as long as 
the time scale r is sufficiently large. We found a strong 
dependence between the minimum time scale Tmin and 
Neptune's inclination. Thus, an external constraint on 
the speed limit of Neptune migration may also constrain 
its inclination. For instance , the u pper boundary obtained 
bv lMurrav-Clav fc Chiallgl (|2005f ) is t < 10 Myr. In that 
case, our results show that under the hypothesis of Sec- 
tion[31 the initial inclination of Neptune's orbit must have 
been larger than 7 deg. On the other side, in all our 
studied cases, the minimum time scale must be at least 
greater than 7 Myr, whereas iMinton fc Malhotral (|2009f ) 
found T < 0.5 Myr. This contradiction may be raised if 
one considers different evolution laws for the semi-major 
axes, eccentricities, and/or inclinations. In that scope, we 
have given in Section [5] an algorithm to know if Saturn 
can be tilted for any migration law. 
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